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SUMMARY 

The Wentzel, Kramers, Brillouin method was used to determine the ionization 
lifetimes of hydrogen, lithium, sodium, rubidium, cesium, and xenon in an applied 
electric field. These lifetimes were used in the analysis of a theoretical 
plane-diode engine consisting of a grid emitter that field-ionizes the impinging 
propellant. Results show that monatomic gaseous propellants are not suitable for 
field- ion emission engines of this design because of the excessively high voltage 
requirements for ionization. 


INTRODUCTION 

The power expended in the Ionization of a propellant Is a primary source of 
Inefficiency in electrostatic rocket engines that employ contact Ionization or 
electron bombardment. An important improvement in engine performance could be 
effected if Ionization could take place at ambient temperatures and nearly 100 
percent propellant utilization efficiency. The idea of field-emission ionization 
arises naturally from these considerations. It is the purpose of this paper to 
examine the feasibility of applying this effect to electrostatic rocket engines. 

The lifetimes of hydrogen, lithium, sodium, rubidium, cesium, and xenon in a 
high electric field are calculated by using the Wentzel, Kramers, Brillouin 
(W.K.B.) approximation. Image potentials of the ion and electron caused by the 
anode surface were neglected. Omission of these image forces can be justified if 
ionization takes place relatively far from the surface. A coulomb attraction be- 
tween the outer electron and the remaining atom was assumed for all the elements; 
this is correct except for xenon. The resulting error for xenon is minor since 
the prominent terms in the W.K.B. approximation axe the field strength and the 
ionization potential. 

The Bohr theory is used to calculate the frequency of the electron in its 
orbit, that is, the number of co ll isions per second that the electron makes with 
the potential barrier created by the superposition of the coulomb field and the 
external electric field. 

The effective capture diameter of an infinitely long, charged, cylindrical 
wire is computed for a co ll imated beam of neutral atoms of given polarizability. 
The Ionization probability is calculated for these induced dipoles, whose force 
of attraction is into the high field region, and yields the relation between 
field strength and cylinder radius necessary for ionization. 
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A theoretical engine is analyzed in which the emitting surface is a plane 
grating of fine vires. It is shown that the space-charge-limited current is 
equivalent to that of a plane diode in which the grid acts as the emitting plane. 


THEORETICAL BACKGROUND 

Two types of phenomena are associated with the introduction of an atom into 
an electric field. First, the degeneracy of the energy levels (with the same 
principal quantum number n hut with different orbital quantum numbers Z) can 
be removed, giving rise to the Stark effect. Second, If the fields are high 
enough, the phenomenon of autoionization, or electron tunneling through the po- 
tential barrier, will take place within a finite length of time. The quantum 
mechanical formulation of these related effects will now be very briefly re- 
viewed. 

The Schrtfdinger equation for a hydrogen-like atom (neglecting spin-orbit in- 
teraction) in a constant external field Is 

- W ** + (- 4 ^ - ^ (1) 

(Symbols are defined in appendix A, and physical constants are presented in 
table I. ) Written in parabolic coordinates, 
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equation (l) separates into: 
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Here ty(|,i},cp) = L(|)H(r|)$(cp), and Kp, Kg, K3 are separation constants. 

Notice that equation (3b) contains terms that go to -00 for infinite values 
of the coordinate Reference 1 shows that such an equation gives no 
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quadratic ally integrable solutions and contains no stable stationary states; that 
is, the effect of the electric field is to give aperiodic wave functions, and 
hence a definite probability of finding the electron at large distances from the 
nucleus. This seems contrary to what is observed in the Stark effect, namely, 
definite spectral lines indicating periodic wave functions that give a vanishing 
probability of finding the electron at infinity. In reference 2, however, the 
lifetime of this unstable state is calculated for the case of hydrogen and has a 

value of 1010 10 seconds for a field of 1 volt per centimeter. The lifetimes are 
still extremely long for the fields encountered in the Stark effect, and an ob- 
server need never be aware of the autoionization phenomena. Spectral lines due 
to transitions with higher principal quantum numbers have been seen to disappear 
under the influence of an electric field. 


EQUIVALENT ONE-DIMENSIONAL CASE 


Because of the obvious difficulty in finding exact solutions to equation 
(l), simplifying assumptions are invariably made. A good approximation (see 
refs. 3 and 4) of the ionization time is obtained by solving the one-dimensional 
problem, namely 


ft 2 d 2 ^ + 
2m e dz 2 



e^z W = E\[r - -e&Jr 


(4) 


and applying the W.K.B. approximation (see ref. 5). The situation is depicted in 
the following potential diagram, which shows the external field superimposed on 
the coulomb field: 



(a) 
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The penetration through the potential harrier is then giyen by 
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The lifetime is found by taking the reciprocal of the penetration probabil- 
ity D multiplied by the frequency v of the electron striking the barrier, 
that i s 


T = 


vD 


( 6 ) 


The limits on the integral give the end points of the barrier as shown in 
sketch (a). Appendix B gives the solution of the integral and the ionization 
time, in mks units, as 
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K(|,k) is the complete elliptic integral of the first kind, E k) is the 

plete elliptic integral of the second kind, and g is the ionization potential 
in volts. 
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It is at 

frt 

this field strength that the barrier disappears and t = v"-*-. 
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A maximum field strength can be defined by letting = 

^ v rmxx 0 


The orbital frequency v is obtained by recourse to Bohr's theory, which 
gives for hydrogen and the alkali metals 


m e e ‘ 


v = 


n 3 2rth 3 (4it e 0 )‘ 


(8) 


Table I gives the values of the principal quantum number n and the ion- 
ization potential £ for the elements considered. Figure 1 shows the logarithm 
of the lifetimes plotted against field strength of these elements. 
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Logarithm of ionization lifetime, log 


TABLE I. - FUNDAMENTAL AND DERIVED CONSTANTS 
(a) Miscellaneous properties of elements discussed 


Element 

Chem- 

ical 

symbol 

Prin- 

cipal 

quan- 

tum 

num- 

ber 

n 

Hydrogen* 5 

H 

1 

Lithium 

Li 

2 

Sodium 

Na 

3 

Rubidium 

Kb 

5 

Cesium 

Cs 

6 

Xenon c 

Xe 

- 


ability, point of dipole at 
a, (at 760 T. and 

(f)(m2) mm Hg), _ . ,• 


13.595 1. 673x10 ' 

5.390 1.151xl0“ 26 1.3x10" 39 1609 

5.138 3.82>a0" 26 3. 0x10“ 39 1153 

4.176 1.419X10" 25 5. 6x10“ 39 973 

3.893 2.205X10' 25 4.6xl0“ 39 943 

12,127 2.18X10" 25 4. 6x10* 40 166 


(b) Physical constants 


Avogadro’s number, Nq, 
atoms /kg atomic vrt 

6.0251X10 26 

Electronic charge, e, coulombs 

1.602X10- 19 

Electronic mass, kg 

9. 11x10“ 31 

Planck's constant, h, (j)(sec) 

6.6237x10-34 

Boltzmann's constant, k, j/° K 

1. 38x10-23 

Vacuum dielectric constant, €q, f/m 

107/4kc 2 


^Monatomic hydrogen. 

c Orbital frequency assumed to he 10 16 sec for xenon. 








SINGLE-WIRE IONIZER 


There are two problems connected with the construction of an efficient ion- 
izer using field-ion emission. First, a very high field must be established by 
some combination of electrode geometry and applied voltage. Since there are 
usually other restrictions that limit the voltage, the high field must then be 
produced by employing electrodes of minute proportions. This Is somewhat un- 
fortunate because the strong fields now exist only in the rather restricted re- 
gion of the immediate vicinity of the electrode, This presents the second prob- 
lem. The incoming gas must be made to remain in this "ionization zone" for a 
sufficient time to allow ionization to take place. The situation is remedied 
somewhat because of the dipole moment induced In the particles by the field and 
their resultant attraction into the strong field region. Since this energy of 
attraction is proportional to the square of the field strength, It Is obvious 
that a point electrode would produce a shorter range force (~l/r^) than a wire 
electrode (~l/r 3 ) and that the former would give a much smaller capture cross 
section for ionization. For this reason this section considers only the wire 
emitter. 


r 


a 


The electric field surrounding an infinitely long cylindrical wire of radius 
can be written 




Va*a 

r 


( 9 ) 


where is the field strength at the wire surface. 

The potential energy U of an Induced dipole In such a field is 


U = 



-a 


r2v.2 


r a" a 
2r 2 


( 10 ) 


Here a Is the polarizability; values are given in table I. In general, the 
differential equation defining the orbit of a particle In a central field Is 
given by 


de 


dr 


r 2mpE 


2mpU 

~w 



( 11 ) 


where l is the angular momentum of the particle; E, the particle energy, Is 
equal to Its thermal energy of translation and is related to the angular momentum 

*>y 


l = s-^ 2mpE 


( 12 ) 


where s is the Impact parameter. 
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Substituting equations (10) and (12) into (ll) results in 



Integration then yields the equation of the orbit: 



The minimum value of r 



( 13 ) 


(14) 


If r ^-j n < r a , then a collision with the wire occurs. The resulting relation be- 
tween the maximum energy and the impact parameter for a collision is 


\ a ^1 r a 


-max ? ? 

a 


Assume that the wire Is placed In a collimated, neutral particle beam of 
flux density whose energy distribution Is Maxwellian, The differential 

flux density dJjj with an energy between E and E + dE is 




The number of particles per second per unit wire length having an impact 
parameter between s and s + ds and with an energy lying between E and 
E + dE is then 2dJ^ ds, the factor 2 occurs because both the top and the 
bottom of the wire are being considered. The total number of collisions is then 
found by Integrating equation (16) over E up to E max and over s > r a and 
then adding to this the number of particles having an impact parameter less than 
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Carrying out the integration over E leaves 



equation (18) becomes 



= 2%R a 

vhere R a is the effective capture 
shovn in figure 2. 



( 20 ) 

against e is 


radius. A plot of R a / r a 



Figure 2. - Ratio of effective collision radius to wire 
anode radius as function of field energy of induced 
dipole for an infinitely long wire anode. 
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The probability that an atom will be ionized on passing through a region of 
high field strength is given by 


P = 1 - 



( 21 ) 


A conservative estimate of the ionization efficiency of a beam of particles 
can be made by considering only the fast-moving ones in a direct collision with 
the wire. These particles will have the smallest residence time and hence will 
define a minimum ionization probability. 

For an atom whose impact parameter is zero (no angular momentum) the veloc- 
ity through the field is 

= J— (E- U) (22 

dt fnie 1 ' 


where U is the potential energy of an induced dipole given in equation (10) and 
E Is the initial thermal energy of the atom. Therefore, 


dr 

dt 



(23) 


The field surrounding the wire is given in equation (9), and differentiating 
gives the relation between dr and 


dr 


*a r a fg 

? 2 


(24) 


In terms of P the 


transit time becomes 

•^a r a 






(25) 


Substituting equation (25) into equation (21) gives the ionization proba- 
bility entirely in terms of 



Setting the integral equal to unity gives an ionization probability of 63 
percent. The wire radius necessary to achieve this probability is then deter- 
mined by integrating over the field strength from zero to 
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The integral in equation (27) was computed numerically for E = 4kT at the 
boiling point of the gas at atmospheric pressure; the results are given in fig- 
ure 3. 1 Over 98 percent of the particles have energies less than 4kT, according 
to the Boltzmann distribution, and of these at least 63 percent will be ionized 
in a direct collision if the relation in equation (27) is realized. For indirect 
collisions the particles remain in the field a longer time and, therefore, will 
also be ionized. 


GRID IONIZER 

The potential due to the uniform grid of wires (shown in sketch (b)) with a 



^Hydrogen is not included in this calculation because its lifetime as plot- 
ted in figure l(a) is for the monatomic gas rather than for its natural molecular 
state. 
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Anode radius, r a , microns 




( 28 ) 


grid spacing a along the x-axis is given by the real part of the function 
(ref. 5) 

^ f ) 

z = x + iy 

The potential is, then, 

V = V a - K In 4 ysin 2 ^ x + sinh 2 ^ y^ ( 29 ) 

The field strength is also obtained from equation (28) and is 

?= ~ = ~ K cot -2 z* (30) 

dz a a 


¥ = -K In 2 


( 2 


The asterisk denotes the complex conjugate. The constants K and V a are 
evaluated in terms of the field strength on the surface of a grid wire and the 
arbitrary zero of potential. The field strength for small z is 


.r= 


a 




(31) 


At r a the magnitude of the field strength is thus, K = J^r a . The 

zero of potential will be taken at a distance b from the grid, a distance equal 
to the accelerator spacing. Since b » a, 

sinh 2 - b = i e 2itb/a » sln 2 £ x 
a 4 a 


and therefore 


V a = ^ r a ~ b (32) 

In order to retain the previous single-wire analysis in its application to 
the grid surface, it will be necessary to space the grid wires so that the local 
electric field about any individual wire behaves as a l/r field at least up to 
r = R a , the effective capture radius. The circular sine and the hyperbolic 
sine functions are equal to each other for small arguments and can be approxi- 
mated by the value of the argument. Thus, for 0 = 0. 3, sin 0 = 0. 2955, and 
sinh 0=0. 3045. 


The equipotentials given in equation (30) are then early circles for 

2 2 

sin2 | x + sinh 2 “ y = (|j) U 2 + Y 2 ) = r 2 < (0- 3) 2 3 r < 0. 3 (33) 
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When r = R a , equation (33) will be made to hold and hence the grid spacing 
will be fixed in terms of the wire radius. The value of from equation 

(34) is, therefore, 

a min= 10 * 5R a < 34 ) 



Figure 4 shows a plot of the grid voltage, given by equation (32), as a 
function of wire radius. The condition in equation (34), as well as the data 
given in figures 2 and 3, was used in calculating these potentials. The accel- 
erator spacing b was chosen to be 2 millimeters, which is probably a practical 
minimum for the plane-diode engine. 

Equation (34) is also used to yield the angle that a collimated beam of 
neutrals should make with the normal to the grid plane so that all the particles 
in the beam can pass within a capture area defined by R a . The relation between 
the effective capture radius and the angle of incidence of the neutral beam 

2R 

0 > cos -1 — - £= 79° (35) 

a 

can be derived from sketch (c). 

A problem peculiar to the ionization of a gas by field emission arises when 
it is desired to calculate the maximum current density attainable. For space- 
charged-limited currents the boundary condition is such that the field strength 
vanishes on the emitting surface. This condition can never be applied to the 
isolated wire or point emitter in field emission, for ionization quickly ceases 
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when the field drops helow a certain level. The situation is somewhat different, 
however, for a grid of wires. Equation (30) shows that the potential of the grid 
surface very rapidly approaches that of a plane. It is estimated in appendix C 
that the saturation of this plane diode results in only a 3.7-percent saturation 
of the individual grid wires. The local fields about the individual wires should 
hardly be affected. 


Grid normal 



(c) 


CONCLUDING REMARKS 

The voltages necessary for the successful operation of an engine of this 
design that employs the field-ion emission technique seem to be beyond present 
capabilities, at least for monatomic gaseous propellants. Besides the problem of 
attaining high voltages, a related problem presents itself, namely, the high spe- 
cific impulse that results from these low-mass particles being accelerated to an 
excessive velocity. Theoretically this problem could be alleviated by the use of 
an accel-decel system. Heavier particles such as molecules or colloids would 
also give a lower specific impulse and would therefore be more adaptable to a 
field-emission engine though the former would introduce fragmentation problems 
in the high field. The greater polarizability and the lower thermal velocities 
of the colloids will result in somewhat larger cross sections for ionization. 
Perhaps further study will show that field emission is a workable means of charg- 
ing these particles. 

Despite the pessimistic results for the particular elements studied in this 
report, it is felt that field-emission ionization is still an attractive method 
of producing ions. The fact remains that ions can be produced without excessive 
difficulty by using needle-shaped anodes and only moderate voltages. The basic 
problem in electrostatic engines will be containment, that is, keeping the neu- 
tral particles in the vicinity of the high field for a long enough time. 


Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, September 27, 1962 
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APPENDIX A 


SYMBOLS 


A 

a 

b 

C ± 

D 

E 


e 

ft 


N 


J 

K 

Ki,k 2 ,k 3 

k 


L(6),H(n),®(q>) 

l 


me 

mp 


atomic weight per ionic charge 
potential turning pointy grid spacing 
potential turning point, accelerating distance 
function defined in eqs. (B2) and (B3) 
penetration probability 

initial thermal energy of atom; complete elliptic integral of 
second kind 

electronic charge 

ionization potential 

applied electric field 

electric field at cylindrical anode surface 
Planck’s constant divided by 2 jt 

space-charge-limited current per unit length for diverging 
cylindrical flow 

flux density of neutral beam 

plane-diode current per unit length of grid emitter 
current density 

a constant, complete elliptic integral of first kind 
separation constants (see eq. (3)) 

Boltzmann constant 

functions defined in eqs. (3a), (3b), and (3c) 
angular momentum and orbital quantum number 
electronic mass 
particle mass 


N 


number of particle collisions with emitter wire per second 



n principal quantum number 

P probability of ionization 

R a effective capture radius for ionization 

r spherical coordinate; range force 

r a vire anode radius 

s impact parameter 

T temperature 

t time 

U potential energy 

u dummy variable 

V electrostatic potential 

V a potential at anode 

W complex potential 

Z atomic number 

x,y, z Cartesian coordinates 

y function defined in eq. (19) 

z x + iy 

a polarizability 

p Langmuir function for cylinders 

e field energy of dipole in units of kT 

€q dielectric constant of a vacuum 

field energy of dipole in units of kT at boiling point and lv/A 

0 angle of incidence of neutral beam 

9 azimuthal angle 

[i reduced mass 

v orbital frequency 
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l;T],cp parabolic coordinates 

*t ionization lifetime 

y wave function, solution of SchrtSdinger * s equation 

Subscripts: 

max maximum 

min minimum 



APPENDIX B 


SOLUTION OF TRANSMISSION PROBABILITY 

According to the W.K.B. approximation, the lifetime of an atom is obtained 
by combining equations (5) and (6) and carrying out the integration: 



The roots of the integrand are: 



With Z l/2 = y and equations (B2) and (B3) substituted, equation (Bl) takes 
the form 



This integral can be -written in terms of elliptic functions (ref. 6, p. 57, 
formula 218.11) 



/ Jt 1 

where sn and cn are the elliptic sine and cosine f- unctions and = Khp,k 
is the complete elliptic integral of the first kind, in which 
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Integrating equation (B5) leaves 

(b - a) 1 L _ j c 2 )E(u) - 2(1 - k 2 )u - k 2 sn u cn u dn u] K 

-y/b 3k 4 L JO 

= ib 1 / 2 [(a + b )E (|,k) 

In terms of C + , 

■fo ~ [(2 - k 2 )E(u) - 2(1 - k 2 )u - k 2 sn u cn u dn uj 

V b 3k 4 L J0 

i e 3 / 2 


2aK(|,k)] 


3^/2 & 


ci /2 [ E (|,k) - C_K(f,k)] 


and the final form of the lifetime is 


T 


exp 



3h 


£ 3/2 




(B7) 


(B8) 


(B9 ) 
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APPENDIX C 


SPACE- CHARGE-LIMITED CURRENT OF GRID EMITTER 

It is desired to calculate the limiting current "between the two cylinders 
r a and R a and find its fraction of the limiting plane-diode current per wire. 
The limited current density for a plane diode is given by Child’s law (ref. 7): 

J = 5.467xlO- 8 A _1 / 2 v|/ 2 b- 2 amp/m 2 (Cl) 

The current per wire per unit length of grid is then 

Jp = 2Rgj ( C2 ) 

Substituting equation (Cl) into (C2) yields 

Jp = 1. 09 3x10" 7 A" 1 / 2 V a / 2 b _2 R a (C3) 

The space-charge-limited current per unit length, for diverging cylindrical 
flow is given by the Langmuir relation (ref. 7): 

J c = 3.43xl0- 7 A _1 / 2 (AV) 3 / 2 R a 1 |3~ 2 (C4) 

The ratio of these two currents measures the degree of current limitation 
of the grid wires in terms of a space-charge-limited diode: 

Jp V|/ 2 p 2 R 2 

J c “ °' 319 AV 3 / 2 b 2 (C5) 


The term AV is the difference in potential between r a and and 

AV = In ^ ( C6) 

x a 

The potential V a is defined in equation (32). The worst case occurs when (3 
and R a are large and b is small; then Jp/J c is at its largest. 

With b = 2X10" 3 meter, r a = 10" 5 meter, Ra/r a = 10, and p si, it is 
found that 



0.037 


(C7 ) 


The amount of current needed to saturate the plane diode is at most only 3. 7 per- 
cent of the current needed to saturate the individual grid wires. 
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